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distortion amplitudes in x and y direc- 
tions of Fig. 1 » » « 

distortion amplitudes in 1,J,K Direc- 
tions of Fig. 1 

combination of terms defined by (32) 
coefficients in (58) 

constant in a convergence producing term 
of 

coefficients in description of subsonic 
vortical flow 'n (19) 
chord 

mean circumferential distance between 
blaoes. Fig. 1 

combination of terms defined by (48) 
and (68) respectively 
arbitrary function 
blade height 
integrand of (54) 

unit orthogonal vectors in upstream Dis- 
tortion coorcinate system 
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♦Aerospace Research Engineer. 

Ifhe velocity ano pressure perturbations 
downstream of the shock are the Differences between 
actual and mean downstream values, however they are 
non-oimiensional iaed by the upstream quantities 
anc as discussed on pg. 574 of [2]. 
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branch points in a space 
combination of terms Defined by (23) 
combination of terms defined by (24) 
wave numbers = MuUy/B^ and 
Mc“'c/Bfc respectively 
spanwise wave number » »q/h 
mean circumferential distance along blaoe 
row 
lift 

Mach number 
moment 

integral number of distortion transverse 
wave lengths in distance s^ along % 
amplitudes of Pn.Pa 

pressure perturbation non-dimensionalized 
by ru^u 

Fourier incices in (] ) 
combination of terms used in (lut) 
combination of terms used in (30) 
solidity « chord/blace gap « 

1 /cimersionless d 

interblade Distance along y in Fig. 1 
(sin x)/S as in Fig. 1 
time non-dimensionai izeo by multiplying 
by <*u/<^ 

combinations of terms defined by (33) 
ano (34) 
mean velocity 
amplitudes of u,v,w 
total. vi«'tical. and acoustic velocity 
perturbaiions non-d imensional ized 
by 

hlade coordinates (Fig. 1) non- 
dimens ion? 1 izeo by blade chord 
X + sY and x + s^ - 1 respectively 
orthogonal coordinates in a fixed 
frame upstream of cascade 
Fourier transform variable 
combination of terms Defined by (69) 

'2 \1/^ / 2\1/^ 

^M^ - ij and (] - M^j respectively 
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comMnation of terms deti.ned b> (50) 
and (67) respectively 
ratio pf specific heats of an ideal gas 



■ •‘ri * respectively 

combination of terms defined by (41) and 
(56) respectively 

imaginary part of integration contour 
for *u 

Kronecker delta; - 0 if i ^ j; • 1 if 
i • j 

small damping terms in yu 4no >d 
respectively 

combination of terms defined by (6?) 

arbitrary variable illustrating perio- 
dicity relation 

variables defined by (43) and 
(60) respectively 

combinations of terms defined by 
(40) and (55) 

combinations of terms defined 


by (1?) and (53) 

angle between upstream and blade coor- 
dinate systems of Fig. 1 
direction of transport of upstream distor- 
tion with respect to axis of 
turbomachine 

combination of variables defined by (4o) 
density 

interblace phase angle defined by (3) 
amplituce of « 
velocity potent ial 
stagger angle « v - u 

(“u * ^g) i"d *^g) respectively 






vl/? , 


I and Ik 

■0 • 


respectively 

arbitrary function in (13) 
frequency . Prp^f^/L in dimensional 
form 

reduced frequenc ies vC|^^, and 
i"C/%o respectively 


SUBSCRIPTS: 
b blade 

d.u downstream and upstream of shock 

respect ively 
n integer 

p.q Fourier indices describinc upstreani 

distort ion 

s shock ^ ^ 

1,?..’ correspond to T.U.k and tc x.y.a 

direct tors in upstream distortion and in 
blade coordinate systems of Fig. 1 
re'pect ively 

SUPERSCRIPTS: 

d.u downstroam and upstream of shock 

respect ive 1y 

upper anc lower half planes of a spare 
' implies both vortical ano acoustic 

propert ies 
vector 

» unit vector 


INTRODUCTION 

Associated with high performance axial fans and 
compressor^ of aircraft engines are slender blades 
with high tip Mach numbers fl]* Such blade rows may 
be vulnerable to distortion of the oncoming airstream 
such as atmospheric gusts, wakes from upstream struts 
and guide vanes, and maldistributions in air inlets 
and ducts; all of which tend to force the blades into 
vibratory motions. Of concern is the influence of 
strong shocks occurring in blade passages. When such 
shocks couple to the distortion pattern they can have 
an oscillatory influence on the blare pressure dis- 
tributions. forces, ano moments. The influence of 
the blade row on the distortion is also of importance 
since it passes on to influence other components of 
the propulsion system. 

The recent linearised theory of Goldstein, 

Braun, and Aoamcsyk [?]. which analyses unsteady flow 
in supersonic cascades with strong in-passage siiocks, 
is applied to the Distortion problem in the investi- 
gation reported herein. The flutter boundary condi- 
tions of [?] are however replaced with those of the 
t hree-dimensional distortion of Goldstein [3]. This 
disturbance profile (Fig. 1) has components parallel 
to. and in two cirections transverse to. a mean flow 
velocity ♦. Relative to the cascade the mean ve- 
locity, ‘Ky. is supersonic in this study. The com- 
ponent of parallel tc the axis of the turbo- 

machine is, however, subsonic. The Mach waves from 
the leading edges of the blades then extend forward 
of the edge of each successive blade permitting in- 
teraction between blades through the supersonic, as 
well as the subsonic, flow media. This is usually 
the case of interest in advance type fans and com- 
pressors. 

The cascade (Fin. 1) is envisioned as an “un- 
rolled annulus" which translates with a blade veloc- 
ity ♦ Integral multiples of distortion wave 
length must however equal a mean circumferential dis- 
tance along the face of the blade row. Spanwise flow 
distributions are considered at limited blade aspect 
ratios although centriiugal and Coriolis forces are 
neglected. Within the confines of lineari7ed theory 
blade thickness, camber, and angle of attack make no 
contribution to the unsteady forces (Chap. 3 of 
[3]). Such time independent effects ran therefore be 
superposed on the results of this study in which the 
bladfx take the form of flat parallel plate'. 

Separat ion of the flow into two distinct (sub- 
and supersonic) regions permits solution of the 
potential part of the flow field in each one by 
adaptation of the Wiener hopf technique [?]. In the 
first region a supersonic flow extends from upstream 
infinity to a cascade which extends to downstream 
infinity. In the second region a subsonic flow is 
guided by parallel blade passages from upstream 
infinity, past the trailing edges, and on to 
downstream infinity. The supersonic solution is 
indeper.dent of downstream influence and the subsonic 
solution is kept sufficiently arbitrary to satisfy 
shock boundary condit ions. Matching of solutions at 
the shock interface involves inversion of a large 
nearly diagonal matrix. This is easily done 
numerical ly. 

In the present effort the shock is assumed to be 
located an arbitrarily small distance inside the 
leading edge of the blade so that it doesn't disturb 
the supersonic flow on the upper surface in this 
region. 
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FORMULATION 


(4a) 


B> the splitting theorem (chap. 5 of [3)) any 
lineari/eO description of an inviscld compressible 
fluid car be decomposed into two non-interacting mo- 
tions; one of which is vortical and the other acous- 
tic in nature. The vortical motion is solenoidal, 
and its substanfial derivative is epual to aero. 
Shear layers, for example, are vortical in nature, 
and are an effective means of forming flo;< distor- 
tions. The acoustic motion, on the other hand, is 
irrotational and has a velocity potential which sat- 
isfies the conveefeo wave eguation. 

Vortic a l fat ion Upstream of Shock v 

The upstream flow distortion of interest herein 
is vortical. Its velocity (fig. 1) is specified by 
the double Fourier series (char. S of [3]) 

E ^<‘*^po ^ Ha’ 

P.d 


a • -.fl^pqSin g * .‘jfpqCos g, 

A . j^pqCOS g .ifpqSin g, (4b) 

and, by the solenoidal condition on v„, 

t • -?ipl .if/OL cos V • - ■ - - f - (4 cos g ♦ A sin g). 

fe L cos V 


(4c) 

Since final results can be superposed, only one set 
of Fourier components will be carried through the 
analysis; thus the p and g indices will normally 
be omitted from the flow variables. 

Acrusl ic Mflirn Upstream of Shock 

The total motion. By ■ vj, « w^, will satisfy the 
irviscid tangency condition a1 the blace surface if 
the norm.al component of ti.e acoustic (potential) mo- 
t icr sat isf ies 


„ i? *pY, II cos V 

♦ KCp^s^n (-aVyhlTe . (1) 

For each set of Fourier inoices the distortion ampli- 
tudes ere j/pq, Vnq. and ^pp ter motion in the 
three orthogonal direct lors i. *, and K. As in [31 
there is an integral number, p, oi dis.ortion wave 
lergths parallel to the face of the blace row where 
L is a niean circumferential distance. The spanwise 
wave length is the hub-t o-shrour distance, ti, divided 
by Fourier index a. Uistence and velocity are non- 
cinif r«ional iaec by blaee ciioro arc near upstrean ve- 
locity. of fig. 1. respectively. 

Transf ormat ior of (1) to blare coordinates 
(x.y,?) movirg with velocity i/j, relative to 1i 
Gives ti>e vortical motion a harnonc tine cepenoence. 
In this system 



111' 

?> 


io (x-*s^-*ns cot g-t ) 

• ’ cos (fcq7)e 


( 5 ) 


tor -sf < X - ns'^ < 0, 0 < a < h, and n • 0. *1, 

•••. By inspection, the solution to the lineariaed 
corveetpd wave eouation, 

( 7 ? - . 0 ( 6 ) 
which sat isfys (*') is of the form 
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where time ha' been non-dimensioraliaeo by tfg/r 
arc whfrf 
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Note that the a comprrents oi tiott tgs. (i) and (8) 
(ab) ifdivicually satisfy the flow tangency condition at 
Fcth the hub arr 'hroud. 


arc spanwise wave numter, fcq, is eaual to *o/h. 

The exponent m (a'b) «fts ti'f irtert.late phase arole, 

r. ^4 


hressure. pj,. 's assoc iateo with the aioustic 
motion, satisfies the wave eouation, and can also be 
expre"fd in ttie torm 


n • Wilis'*^ • S cot g). (3) 

Ihc X and t cepencence in (aa) arc (ab) together 
sati'ty the corcitior that the suh'tartial rcriva- where 

t ivt , ("Vy/tit, is eoual to aero. The cotrpor»r-ts ct 
the cistortior amplitudes in the x,y,a rirectiors 
arc row 


-ig. t 

pj^ ' -a f “ cos (fcqa)h^j(x,y) 



(‘•a) 


(yt) 



ORIGINW. * 

OF POOR QUAtlTY 


The prpceCiire for spluticp of ♦u(>'.y) by use of 
the Wloner-Hopf technique Is given in Appendix 1. 

The results for 0 < y < s ere 
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for X < - a^(s - y) (10a) 
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Siihst i tuf ion of no) into (8) arO (9) detertrincs the 
velocity of the acoustic motion and the pressure in 
the supersonic regicn ahead of the shock. 

Vort i c a1 Motion [downstream of Shock 

Tn this region Tet the velocity, V(|, take the 

form 


which satisfies the spanwise wall tangency condition 
and has t substantial derivative equal to zero, lo 
satisfy periodicity let 


iw.(x«ns^) 

!> S(y ♦ ns) 


j_ iuj* 

e^"" e S(y) 


(14) 


for blade number n • 0, *1, *i, The sclenoi- 

dal condition will be invoked later, completing the 
splitting theorem requirements for V(j to be vor- 
tical. ^ 

It is anticipated that a part, Oj], of vp can 
be found which is both vortical and acoustic, as in 
to. (5.1?) of [3]. Coirtiining V() with the poten- 
tial flow, W(|, may then simplify determination of 
the remaining part of V(j. Such a procedure^is 
used to advantage in [}■]. Associated with v^ is 
a potential *() such that v>,) • 9o, and sinc^ 

V • vp . 0 it follows that »''ert « 0. Since vp 
is vortical, fp makes no contribution to the pressure 
distribution. Let fp take the form of (13) 

i(“rtX-%,t) 

. a cos (fcjjZ)e f(y) (16) 

For fo to satisfy Laplace's equation, f must 
satisfy 


■("d " ^^o)^ F" • 0. (16) 

A solution of (16) which enables to satisfy 
periodicity as well as to cancel the normal component 
of Vji at the blade surfaces is. 

Op (ns) r 

’ ■ T vrn - h t t; t ) f’o<> - 

- cosh rT^(y - ns - s)t| (17) 


where ifj. * fc' . The subscripts l.?,3 on o 
idrrtify the components of O in the x. y, and 7 
directions respectively. 

Substitution ot (13). (16). and (17) into thr 
icjcrt if V 


Vj. i (Vp - v»,J ♦ (18) 

shovs that the first term (V(i - vjp) makes no contri- 
bution to the Mace tiouncary cordinons (since its 
y component vanishes at y » 0 and s) although it 
does contribute to the matching of flows at ttie shock 
boundary. With sufficient generality the strictly 
vortical part i)f (18) car thus l.e sufticirntly des- 
cribed by Fourier sine and cosine series as 


1 (w^X-ki t ) 

v^ • a cos (fcga)e “ S(y) (13) 
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2^ b, _cps(m».v/s) 
m.1 


S bj, 

’ 

-fcntar (fc^7) £) b, „,ccs((tH.y/s)| 

(19) 


'-0-' 


arci the v»o part pt (18) nill he Induced in 
‘'d * ’*d' procedure Is close to that of [?1 

inhich utilizes a stream function for a two-d imen- 
sional vortical motion. 

The scleneidal condition has already been in- 
vokec in Usino it in (19) for Vf, dives 


“'o‘’l .m 
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^ I'-i 
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^c pu stic >Vitipn downstream of Shock 

The acoustic field ooKPstream of the shock is 
composed of two parts. An infinite "duct" solution, 
to the conveetpd wave equation, havinq downstream 
running waves accounts for shock boundary condi- 
tions. A "cascade” solution, on the other hand, 
gives the upstream (or reflected) waves. Each solu- 
tion, as well as the vortical motion of the preceding 
section, individually satisfies the wall tangency 
condition. The sum of the two acoustic solutions 
must satisfy the Kutta condition. In both solutions 
the velocity potential, »(,, again takes the form 

- i w t 

f^(*.,v.a) • a cos ^ •(,(>'. v) (?D 


The procedi're to obtain the comt'ined potential 
♦jj . ♦ jgascade Qiypf, Appendix 11. 
The results for 0 < y < s are 



for C < > < 1 - s^ 


(iif) 




where 
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T, * fcl. i‘ the krcrecker Celti. 

[if t erffinSt lon^ot the cP{ft1cients B|,; ct (??) i<nC 
t'' t (j » 1.?.3) n (19) #rc (?0) tf maer prssihle 
wift the arriticeal irfrrr'aticr supplier hv the 
tpllrwifiq three shock relatiors. 

Relations Across Stiock 

The mear ('low satisfies the Rank ine-Hugoniot 
relation*. Shock curvature arC cisplacement are 
treatec as prrturhat ior< as in [8}. The cownstream 
vcrticitv, loroitufiii.al velocit\. arc pressure arc 
rt later to the upstrean' vcrticitv, lorguuclinal veloc- 
itv. arc. pressure hv U»’ follrvtoQ three eouations: 



(?S) 





(?e>) 


anr 



these eouations Otter trom those ot [irl h> aocition 
ot an upstream vorticity anc a spanwise distribution 
pt pressure arc velocity. 


k*atc hi ng ot Solutiors a t Shock Bo un oary 

fguatior (?b) car easily be cpnFTnec with lt'7) 
to >ielr 



(s’?) 


since the 7 ''ependerce enter* roly through 
cos (fcqa). Crmtiinino (J-R) with (M, (19), anr 
solenoiral ernoition (?C), and utilizing rrthrgonal- 
ity relations between sine terms, permits determina- 
tirn ot bj^p and 83 n ’n term* ot bp p 
and known Quart it les.' the rewn*trean vorfical veloc- 
ity rt (19) can now be expressed as 
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- (-1) exp (iw^s clr u) 
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IT - (u'^S ctr h/t) 


for m > 1 = C 

> 0 for IT e 0. (?0) 

Use of (?). (7), (9). (lOh). (?1), (2?a). arc* (?9j 
for velocities arc pressures in shoc^ '•e’^tlors (25) 
CIO (26) arc utilizing (377). (378). («88). (489) ant) 
(692) ct [9] gives the fcllov.ina eouation for Bp 
ir implicit form. 
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(31) 


where 


i(p»iui,S cot u 
(^) • 

i(r-c s^+ni-) I 

f J!— -jJ 


i'w , A 

*1 ' ~ir ^ 


(32) 


(33) 


♦ J-'u- 

tp * — e 


(f ) • 


(34) 


♦ fef 


The left side of (?1) reduces to (3.33) of [2] as 
feq * 0. The right side of (31) however replaces 
the flutter terms of [2T for this distortion prob- 
lem. Solution of (31) for Ep is by inversion of 
an infinite matrix. The off-riagcnel elements of 
this matrix came from the second term of (31) which 
represents acoustic waves reflected from the back end 
of the cascade. These are rapidly ceceyinq waves 

for (nr/e s)^ ■* ik leJ^ > k^ . so that, 

ir practice. (31) Canute approx iirotnd oy a nearly 
ciagonel matrix. This can be further truncated to a 
reasonably small number of terms, Numerical results 
suggest that as few as ten terms along the oiagonal 
may, at times, give satisfactory results at values of 
recucoc frequency less than one. Once tiie Ep are 
krowi' tiie coefficieit b? p in the Fourier series 
cl (29) car he ottained from shock relation (26) as 


^?,r ' — “d < 




-ir;(l-s') - 

m?0 


?(> ^ ns 


a * rH ^ o'p) 



M K ,\^ 

R 


(36) 


lise of (36) in (29) determines the subsonic 
vortical velocity for 0 e y<^ s and x > 0 wi.ereas 
use of partial derivatives cf (22) in (21) determines 
tt.e ccrresporcinq acoustic velocity components. Use 
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of (10) In (8) <>1ong with (?) and (4) gives the 
acoustic and vortical velocities for 0 £ y < s and 
* < 0. Pressure can be obtained by use~of J^) and 

-Iwyt 

Its subsonic counterpart p. > a e cos (fe.e) 

(1k>^ - »/8*)*j. Extension of equations to other 
blaoes and channels Is singly by use of the perio- 
dicity condition On(x ♦ ns*, ns ♦ y) « e’''*' 
where « can be any of the physical variables 
•u. Vu. Pu- Pri* Vq. and as In [?]. 

RESULTS 

Variables which may be arbitrarily (Independent- 
ly) selected for Input Into the formulation for the 
general case are ratio of specific heats, y. Mach 
number. My, two of the three angles, u, u, and Xt 
cascade parameters S, h, and either phase angle, e, 
or frequency, wy, along with olstorton parameters 
L/p, q, and either amplitudes a and A or j/ 
and Certain combinations of these parameters 

give a very orderly upstream distortion (vortical) 
motion which In turn, through the blade surface 
boundary condition, results In an especially clean 
acoustic motion upstream of the shock. Numerical 
results based on such flows are relatively easy to 
Interpret and should aid In the analysis of moio gen- 
eral cases which are within the capability of the 
theory. 

As In Flo. 2 let v • 0 and -u • x « wave angle 
given by sin*’ (1/M^j). Then -cot • S/S*. Let 

one transverse wave length fit the distance ByS 
along x and let s* « Nb^s for Integer N. Then 
b 2 « 1/N and Ml, . (1 ■* This v*ave pat- 

tern has a negative phase along ay « s which off- 
sets the phase of the cycle along ax > b^s, leav- 
ing r • ?i(N - 1). Using (3) along with d • 1/S • 
s*(] ♦ cot‘u) It follows that wy . ?t(N - 1)S(1 •* 
cot'u)/(l - cot^ u). The lowest N which clear- 
ly permits a subsonic mean axial velocity has a value 
of ?, For N . 2, r . ?T, My . {2l2yl^, anc 
X » -u » tan'l vT. Selecting a solidity of 1.3 
then gives «y - ?0.0 (which Is the passing fre- 
guency) and circumferential wave length L/p • d ■ 
0.769. Amplltuoe ratio A/a Is set egual to 0. 
Amplitude ratio C/a Is then equal to -4.7 1/fen 
by (4c). These results, along with y • 7/6 and 
feq « ./3 will be used as the basic reference con- 
dition for the numerical results herein. 

Influence of Distortion on Blades Pressure Plstrl- 
but Ion 

Fressure distribution along the upper and lower 
surfaces of the blades Is given on Fig. 3(a) at the 
basic crnerec condition. The short Cashed line In 
the supersonic region represents a single airfoil 
solution [10]. Consistent with the reference cascade 
conditions on Fig. ?, the input to the airfoil theory 
Is a transverse distortion having a wave length egual 
to s*/?. The resulting pressure distribution Is 
normallaec to match that of the cascaoe solution at 
the leading edge. The short wave length (high fre- 
ouency) of the Bessel function solution of the air- 
foil Is characterlst 1c cf the cascade solution also.^ 


^Use of a reduced frequency of 20 thus permits 
comparison of two mathematical solutions to the 
linearized converted wave equation (6). It Is 
possible that (6) may net offer a suitably accurate 
description of the motion, however, at such a high 
frequency. See, for example, chapter 1 of [11]. 


In the cascade this short wave length characteristic 
Is carried a short distance downstream of the shock. 
The periodicity relation, ♦(x,y) • e*’*’ *(x ♦ s*, 
y ♦ s). Is used In (2?) to relate pressure on the 
lower surface of blade 1 to the upper surface of 
blade 0. 

A second cascade solution (not shown) for a 
distortion with a phase angle of «/? with respect 
to that of Fig. 3(a) Is obtained from 

r <“u^ 1 

Im[p e /a cos (fej,?)J. This has the same wavelength 
characteristics and^'ln the supersonic region agrees 
with single airfoil theory. Consistent with shock 
curvature and displacement, the pressure perturbation 
amplltuoe Increased for the Imaginary part, but de- 
creased for the real part of the solution across this 
boundary for the reference conditions of Fig. 3(a). 

The solution for Fig. 3(b) Is based on the same 
set of parameters (tor the ordered vortical motion) 
as Fig. 3(a) except o, and therefore by (3), 
are reduced from the basic reference conditions by a 
factor of 1/8. Single airfoil theory Is no longer 
close to these results and has been omitted. Trends 
are quite smooth and, of course, much more graoual 
than those of Fig. 3(a). 

At very low Interblade phase angles, < ^(0.01), 
the pressure profile (not shown) across the channel 
Is gulte uniform and tne pressure amplitude In back 
of the shock on the upper surface approaches that on 
the lower surface near the leading edge. 

Influence of Cascade on Distortion Profile 

Figure 4 uses the same bask reference condi- 
tions as In Fig. 3(a). Envision Fig. 4 as a view of 
the x-y plane of a blade channel. Located In this 
plane by use of the solid sytr.bris a'‘t the leading and 
trailing edges of blades 0 and 1 ano two reference 
Mach waves. Also located at various x stations 
are y distributions of the real part of the solu- 
tion for the three transverse components of velocity 
In up « vp < wp. Far upstream oi the cascade 
the amplitude of the acoustic motion (given by the 
short oashfc line) goes to zero, however the vortical 
motion (long dashed line) is transported undlmlnlshed 
as If frozen In the flow. In the vicinity of the 
cascade the acoustic motion with Its characteristic 
short wave length, cut to diffraction about the lead- 
ing eoges, becomes quite prominent. Once Inside the 
leading edge Mach wave of the nearest blade (blade 0) 
the transverse component of the acoustic velocity, 
wj^), 1s noticeably large only In the vicinity of 
the blai'e Surface. The amplitude of the total mo- 
tion, u^'^1, away from the wall Is then mainly con>- 
prlsed of the relatively long wave length vortical 
motion, At X locations between the trailing 

edoes of blaces 0 and 1 the amplitude of the acoustic 
motion again becomes quite prominent, where diffrac- 
tldr In again encountered to satify the Kutta condi- 
tion. Eventually wJ^J decreases at large x to 
the asymptotic value at the last station, x ■ -. The 
amplitude of the vortical motion continues undlmln- 
Ishec far downstream of the cascade and Is not plot- 
ted at the last station. 

The y clstritjut Ions Of wj*') has much the 
same nature as w^ ''. Although A, ano therefore 
are zero, a longitudinal component cf vortlc- 
Ity velocity, vj®), arises at the shock. Beyono 
this station both acoustic and vortical motions have 
trends with x much like those of wj^) anc 
Recall that part of the vortical motion, 

Vjj has been Included In tg thus slmpllfylnq the 
wall boundary condition. 
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The y olstribut Ion rt «1 sp fo11o*>s 

a pattern somrHhat similar to Ihc Intluence 

of the shock on v\®' Is small out quite pronounreo 
on for the particular set of concutions of 

Fig. 4. 

The a oistrihiit ions of vortical anc* acoustic 
motions are not plotted since, hy {?), (8), M9) and 
(39). they simply vary as cos (kqz) for the trans- 
verse components of velocity and as sin (feg?) for 
the spanxise componrnts in both subsonic arc j uper- 
sonic reqions. 

Influence of (listortion and Cascare Parameters on 
lift and Moment 

Lift, S^, and moment, uf. were obtained hy inte- 
oration of the pressure distributions over the blade 
surfaces. Moments were taken about the center pi the 
blade located at »•!/?- s^. Mots of the imagi- 
nary parts of lifi and rcment v*rsus the correspond- 
ing real parts (chap. 9 of [1?]) can be used to indi- 
cate the influence of various cascade and distort ioi' 
parameters in stability studies. Ihe contribution of 
shock moverent to these terms can be obtained from 
plots as on Fig. 5(a). As in [?], the shocks are 
assumed to be attached to the lower surfaces of the 
blades (fig. 1) arbitrarily near their leading 
edges, the allowed displacemenes , «^, of the shock 
footprints on the upper blare surface, however, cau- 
ses a lift with amplitude equal to (p^, - Pi,)*j 

within the confines of lineariyed theory. Here p;,’ 
and P(< are the steady pressure* before and after 
the shock. This effect 's incluoed in the total lift 
and inoment of Figs. 5(b) and (c). These figures are 
arbitrarily plotted at an intprblade angle of 1/8 
that of the reference set of conditions as in fig. 

3(b). 

The variables h anc a enter the eouations 
for velocity potential (therefore tor lift and mom- 
ent) only as a/h and thus their influence can hr 
pspressro t>> spanwi'r vavc rumter fcg • ec/h. This 
wave nuirher shoulc be most influential on anplituncv 
♦ ,i and when it is greater tfian k„ and k(. 

respective^. This trend is illusfrater in fig. 5 t>> 
the largo distance between tick marks on these curves 
at higi' fcg. Increase of upstream distortion ampli- 
tude ratio, A/a, tencs to increase the magnitude of 
■y and uf at large fcg while preserving the spi- 
ral like (bg dependence, tusps tr\cicr loops charat- 

teriye both the and y curves. The magnitude 

Of y^ is a significant art of y. 

(OMLUl'lNG RfMARkS 

The influence of the third dimension cn this anal- 
ysis is primarily through the a/'muthal wave numPer 
feg. As kg approaches aero the final eguations 
ano their solution reduce to much the same form as 
those of the tliilter analysis of f?1. The main differ- 
ences are in the boundary related term* where blade 
oscillation e>pres*ions of [i] are rppl 'ced bv flow 
distortion esprrss'Ons of the present study. 

Selection of a cuffbinatur of cascade anc distor- 
tion parameters fcr a very orderly irilial vortical 
m.cticr gives in turn, through use of the blade surface 
boundary corcition, an especially orce rly acoust ic mo 
t lor as well. Nuenerical results based on su'h flows 
are relatively ea»y to irterpret and aid in the analy- 
si* cf mere general cases which are within the capatil- 
ity of the theor,. At these special crrcit ions tfie 
supersonic pressure distribution reduces to shapes pre- 
oiefad bv sirule airfoil theory. At high solicity the 
pericdic patterr of the downstream running waves in the 


subsonic regitn of the rascaue Is quite evident In the 
blade pnessune distribution. Intluence ct the re- 
flected, or forward running, waves Is mostly confined 
to within a nelallvely small distance of the trailing 
ee.qes cn the upper surfaces but Is apparent over the 
distance 1 - ?s* < » < 1 - s^ on the lower 
surfaces. ” 

the strong In-passaqe shock has a large effect on 
lift and moment. Numerical results show large influ- 
ences of Inlerbladp phase angle arc therefore reduceu 
freguency, longitudinal as well as transverse clslor- 
t Ion amplitude, and aylmuthal wave number on total as 
well as shock induced lift arc mH'mfnt. this suggests 
that any one of those oistortion paramieters could have 
ronsinerable influence on the forceo vibration effect 
in a « t abi 1 i* y analysis . 


APPENPIX 1 

SOLUtlON FOR AfOUStlC MOtlON UPStREAM OF SMOfK 


Substitution of (7) into (6) gives 



(38) 


A solution by stparatior. of vari'MfS is 
»>p {-ifo - Mi,ku)> f auvu.v]| hheru n is an 

( ? ? ? 

arbitrary constant and > . U - 

Multiplying this result by f(n) e>p [-i(o - M^,ky)s*] 
anc integrating o*rr aimissabli values ot » is alsr 
a solutidn df the linear hrmegeneous Eg. (10), when 
f(o) is ye. to be rntcrmlnec. 

Not, tdllowing the prccedurc Ot hanr anc Frird- 
maf' [At, apply lhi« re«ull tr rarh of r. blades 
using the local cccrdtnatr system ns^, 

Vp . y - ns arr specify that fp(o) must sat isfy 
th» periodicity rone' it lor f?] 

t„(-i) - fr(") (37) 


tor n . 0, *1. *?. ••• . A solution with Jumip con- 
ditions aerdss the r'** blade in proper form is then 


*u.r **n*'n^ 


sgi > I 
"7T- 



i(o-M k 
u u 



((«-M^k 


u 


*'u\ 


l>r|3 


Co 


( 38 ) 

when sgn y • <1 for y < 0. Thi branch points for 
>y are set at *(*y ♦ lt„) wlere . kj •• 
ant fy IS a smi/ll ramplrg terr. Nrw. a* in 17). 
let f • Myiy and Im kj_ . ly. tt'er* along Ifie integra- 
tion fortcur 0* (78) Im (o - M, k, ) . 0, Im >y > C, 
and |f»p{ir(« - S,ky)'sT , ai,Vyjvr|]} » 1. Sumrieg the 


g 


roPtrlbutlon of each blaor to the potential “lelds a 
qroantrlc series which can be enpressed In closed 
fonr as 








-Ko-M k 

y)r “ do 


(39) 


where 


A„(n.>) I 


arc 


7T 


/ - \ 




.A 

VW)\ 


sir a. 


sin t 


u 


tor 0 < y < s 


(40) 


a* I i (p - M k s^ ♦ os^ * ( > s). 

u 7 \ II u u’u / 


(41) 


equations ( 4 ?) and (44) are In the fonn of the upwash 
ana continuity expressions of Ref. [4], They differ 
from those of f?] only by the presence of kg/ib 

1w > 

In ^ and by the distortion pattern, e “ i - 
placiKq the flutter motion on the right side cl 
( 4 ?). The Integrals In both of these equations are 
In the form of inverse Fourier transforms forming a 
two-part boundary value problem which can be solved 
for fq(«) by the Wlener-Hopf technique [6).^ Ut11l2- 
Inq unilateral Fourier transforms f7] over iy, the 
two equations are comblneo and then factored and de- 
composed Into parts which are analytic In respective 
hall planes as In Appendix 0 of f?J, have a comnon 
region of analaticity between them, and vanish at 
•. Ihe solution is coapleteo by use of Llouville's 
theorm giving 



Hfrp ,p Is a damplnq term on the distortion to en- 
sure convergence as Xy • • and can be set egual to 
irre upstream o1 the shock. 

Ihe solution lor ln(o) required that •^,{a) 
could he factored as «„(o.0) • «,*((i) /«„(« ) where 
«*,(o) and «y(o) are analytic and non-yero In the 
upper ( Im a > M,,Cy) and lower ( Im a < M„t„) half 
planes respectively. The factoring procedure of Ap- 
pendix C of [2) yields the results 


Applying t‘.e flow tangency relation (J>) to the sur- 
tace rl the r • blai^e yields 



'•u 1 

JT ’ 77 


(a) • e 




TT (' - “/'’d)’ 


.•nc* 



Me M 

it u 




-1(o-K )x 


U If ' 


Co 


K 

U U 


. W 


U li 


( 4 ?) 


I IT «„ I » s'^ > 0 and vherr 

«. (a.O) I >A. />>! . (43) 

^ ly.C 


Kii piterlial d the irrctat itna 1 field produced up- 
•irean of the cascaof must hi contlruous across 
line? > • ns extftding forward cf lie leaciiig eoges 
(f th» blacis. Applyinq this bouncary condition to 
l?8) ct the /iretb Mao yielcs 



(44) 


- ^<1 


la ) • e 




a y 

u u 


sin (a Y s ) 
u’u ' 
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Equations l^fi) through (60) flUfrr from those of Rf t . 
r?] only by the presence of kg In the oeftntttons 
of end vj. 

Substitution of (45) Into (39) enables determi- 
nation of •i,(*,y) by the method of residues, the 
path of Integration Is closed In the upper half-ptane 
for K < -$y(s - y) and In the lower half-plane for 

* > -luT' - y)* 

ARPENOIt II 

SOLUtlOh FOR Acoustic MOtlON DObNStREAM OF SHOCK 


By (64) and (SB) the lero upwash condition w*ll be 
satisfied If 


/ 


where 


-1(o-M>.)». 

’o n‘“^“o*"’°'' d« . 0 for < 0 


(59) 


An Infinite "duct" sclutlcn to the converted 
wave eouatlon for downstream running waves can be 
ohialned by use of the method of separation of 
variables, this Gives the result 


^duct, 1 


where 


and 


E ’v. 

D e cos 


(61) 

(d) 


n«0 


n* • 

n CO n 


. 1 

n 


[;r./a^s)^ - •'S; ‘ (53) 


(5e) 


ig («,c) I ~*^gi*y (« 

y«0 




-1 


(60) 


A relation speelfylng yero pressure jump across lines 
eatenolrq 1r the downstream' direction from *ne trall- 
Ino edoes of the blades 1$ obtained by substituting 
(61) and (641 through (68) into Pg • ( Iwg - ?/>*) 
^^durt * and utlllying the periodicity 

relation. The resulting eaprtsslon can be reduced tc 


j ’o.n<“’‘' * 


for »j, > 0 


(61) 


By chapter 6 of [3] the cascade potential for up- 
stream running waves can he e»pr*ssed a* 


^cascade, , 1 

•d • 7T 


/ ’o(")*g("> 


y) -<(«*f't,R(,)v 


Solutirn ot Eos. (6®) and (61) aoair constitutes a 
Wlenfr-hrpf prohlen which can be solvtc In much the 
»aae marr*r as In th» sup*rsor’t case yielding 


da 


f (f 
t'.n 




1 ^r^,S^-r*n*J 


(6^') 


(64) 


wharf, in sllohtly different nomi nc la t ure than [3] 
1 


*d • 7T 


r — ^ * ^T- 

sin sm 

L d d J 


(65) 

* Vd»’ * ‘ W''" 

(66) 

‘’JO'" 


(67) 


For tht* linear prof Urn. where •uperpesttU' applies, 
let 


here It i« ass'imer that ran t/e tactored Into 

the form K^(a)/(n(<>) where «f(o) and rfla) ari 
analytic and nor-7*ro m uppar art Irwir half plants 
rrsp< ct Ivt ly. To perform this f at tor 1/at ion suh- 
•titiite (66) into (fO) arc otitair 


(a ' <(') • F..'/. ‘ ^fh 

Odd c c 


(Fd'd'>/(‘ 


'/ «ln a‘ Sin a'] 


•63) 


The numerator of (aj) car of factored by the 
Ueleretras* f ac t or 1/at lon torrxjla [7]. PaOf 40 of 
ff) IS helpful for application of this thiorem to the 
i.'ncminator. Choose the half planes as In fip. 8 of 
(dl but fith (ranch ppifts Of ip at *(*f • Wj) 
vlere . kf - fc(,/rc. Then 




(68) 


11 


«^(o) . 


(» - ♦j,) »1n 


1ob/» 


^nd Ip as 
cosh [a 




(a) $1nh (»_os) mp 


['•(:, - •’)]} 


- CCS (- • Vd*')]0 * -'“o) 

TT 

« I -r — (64) For Ii thf path of Intf^ratlor must bf closed In 
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as a • 


and 




•d<« " *c’ 


Tf 


(l - a/.;) (l . a/«:J 
n.l 


(65) 


Khert 


? 7 7 

tn k‘ orreter than and less tbar •> /»! respec- 

tivel>. Eouaticn* (64) through (69) reduce to those 
of Apprrcl* c of [?1 a! fcq • 0. Substitution of 
(6?) arc (55, through (5P) into (54) gives an irte- 
grel f xpressior: fcr a^*'cadf t-o 

eva1uat(0 by use of the resloue thecrw. The Inte- 
orarc of (54) tar be f»pr*ssec ts 1 • Ij - 
vhere the a ceperrence ef lrrs Ij a* 

cosh r»p»(.v - *)y|»c“^‘c(“) ** 0 "'^ ^'**d|] 


thf upper half place (UHP) lor »c * ® 

Jj] > 0. Ihe path for Ip anst be closed In UHP for 
>0 < s* and LHP for ?o > s'^. Finally, the 
total amplitude, eg, for 0 1 y i s Is obtained by 
comt Inirg aj'^cf vilth a^hscade. 
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(a) Distortion and cascade at basic reference conditions of Fig 
ure 2along with S ■ 1. 3. aiyj • 2ft 0, and • nl\ 

Figure 1 - Pressure distribution on zeroth blade. 
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ibt Distortion and cascade conditions same as Figure tiat. except that o and are reduced to 
2rr/8and 20 0/8, respectively 

Figure 3. - Concluded. 
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conditions same as in Figure 3. 
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Longitudinal distance. x/(l-s’l 
(b) Downstream of shock. 
Figure 4 - Concluded. 
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(a) Shock induced lift. Vj. 

Figure 5. - Influence of azimuthal wave number, d, 
amplitude ratio. A/a. on lift and moment 
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i. and distortion 






